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We study the cold and dense regime in the phase diagram of two-color QCD with heavy quarks 
within a three-dimensional effective theory for Polyakov loops. This theory is derived from two- 
color QCD in a combined strong-coupling and hopping expansion. In particular, we study the 
onset of diquark density as the finite-density transition of the bosonic baryons in the two-color 
world. In contrast to previous studies of heavy dense QCD, our zero-temperature extrapolations are 
consistent with a continuous transition without binding energy. They thus provide evidence that 
the effective theory for heavy quarks is capable of describing the characteristic differences between 
diquark condensation in two-color QCD and the liquid-gas transition of nuclear matter in QCD. 


I. INTRODUCTION 

Tremendeous experimental and theoretical efforts have 
been devoted to the exploration of the QCD phase dia¬ 
gram over many years |T]. Lattice QCD studies were able 
to clarify the nature of the deconfinement phase transi¬ 
tion for vanishing net-baryon density and calculate the 
equation of state to high accuracy mi- A good de¬ 
scription of the lattice data in this regime is nowadays 
being obtaind in Polyakov-loop extended chiral models 
[3]. Based on functional continuum methods for QCD 
such as the Functional Renormalization Group 0E3 or 
Dyson-Schwinger equations [7] these studies can then be 
extended to finite baryon chemical potential [8]. For suf¬ 
ficiently high temperatures one can furthermore use re¬ 
summed perturbation theory to include baryon chemical 
potentials beyond Taylor-series expansions [9]. 

The situation is different, however, in the low temper¬ 
ature and high baryon-density region of the QCD phase 
diagram as relevant for the properties of the inner cores 
of neutron stars. The fermion-sign problem impedes in¬ 
vestigations of this region by standard Monte-Carlo sim¬ 
ulations m- Methods to circumvent this problem are 
under active development: Complex Langevin dynam¬ 
ics has been shown to work for finite-density QCD with 
dynamical quarks at least on small lattices m- Monte- 
Carlo on a Lefshetz thimble has been successfully tested 
for complex scalar fields with chemical potential JT2]. 
Dual lattice representations can help to overcome the sign 
problem for scalar fields and Abelian lattice models with 
worm algorithms hms]- In a similar spirit, graph rep¬ 
resentations in terms of hadron worldlines can be used 
for lattice QCD with unrooted staggered quarks at finite 
density in combination with strong-coupling expansion 
techniques naini- 

To get access to the continuum limit, one can com¬ 
bine strong-coupling and hopping expansions in order to 
derive three-dimensional effective Polyakov loop theories 
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for heavy quarks fT%H20) . It has been shown, for exam¬ 
ple, that such effective theories can reproduce the crit¬ 
ical couplings for deconfinement in the pure SU (2) and 
SU(3) gauge theories, and the critical quark masses lim¬ 
iting the region of the first-order deconfinement transi¬ 
tion in SU( 3) with heavy quarks, within less than 10% 
fl8l IT9l . Studies of the cold and dense regime of the QCD 
phase diagram have produced first results on the nuclear- 
matter transition for QCD with heavy quarks eh i23i. 
In particular, the effective theory yields a binding en¬ 
ergy per nucleon and a nuclear saturation density in the 
right ballpark relative to the heavy nucleon mass, with 
evidence of a first-order liquid-gas transition at very low 
temperatures which are currently still beyond the range 
of applicability of the expansion, however. 

More direct evidence of a first-order transition anal¬ 
ogous to that of nuclear matter has so far only been 
seen in G 2 -QCD at finite density, a QCD-like theory 
with fermionic baryons but without sign problem which 
can therefore be simulated with standard lattice tech¬ 
niques [28l [24] . Another QCD replacement without sign 
problem is two-color QCD with the gauge group SU( 2). 
While it does not have fermionic baryons, it is much 
cheaper to simulate than G 2 -QCD with 7 colors and 14 
gluons corresponding to the fundamental representations 
of the smallest exceptional Lie group G 2 of rank 2 which 
contains the SU( 3) gauge group of QCD as a subgroup. 

The cold and dense regime in the phase diagram of 
two-color QCD with diquarks as bosonic baryons does 
not resemble the properties of nuclear matter but rather 
those of QCD at finite isospin density with charged pion 
condensation [25]. Detailed lattice results for the phase 
diagram of two-color QCD are available [25H2B] and its 
qualitative features are captured by effective chiral mod¬ 
els rather well [29]. They reflect Bose-Einstein conden¬ 
sation of the baryonic diquarks in a second-order phase 
transition which occurs at zero temperature when the 
baryon chemical potential ps reaches the diquark mass, 
and the analog of a BEC-BCS crossover for larger values 
of b as chiral symmerty gets gradually restored inside 
the diquark superfluid. 

In this paper we apply the strong-coupling and hop- 
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ping expansion techniques of Refs. m m to two-color 
QCD, in order to test to what extend these relatively well 
established features of the cold and dense regime of the 
two-color phase diagram are reflected in the effective the¬ 
ory for heavy quarks. Our results from the next-to-next- 
to-leading order in this expansion scheme indicate that 
the corresponding three-dimensional effective Polyakov 
loop theory is capable of describing characteristic differ¬ 
ences between two-color and real QCD at finite ■ One 
can furthermore directly compare our results to lattice 
simulations of the full theory in two-color QCD at finite 
/is to assess where the effective theory eventually breaks 
down when going to smaller quark masses in the future. 

II. THEORETICAL BACKGROUND 
A. Effective Polyakov Loop Theory 

We start with the derivation of the effective action. 
We will only briefly summarize the steps leading to the 
effective action. As most of the steps are very similar, 
we point the reader to a very detailed derivation for the 
the case of SU(3) given in [22]. We start from 3+1 
dimensional two-color QCD with Nf flavors of Wilson 
fermions. After integration over the fermion field, the 
partition function reads 

r N f 

Z= [dUj expb^HdetDy, (1) 

J /=! 


The easiest way to determine the contributions of S\ is 
by using character expansion of exp[— S g ]. In leading or¬ 
der this results in a nearest-neighbor interaction between 
Polyakov loops 

S[ L0) =\(p,N t )Y / L s L y ~, (6) 

(xy) 

A(u(/3), N t ) = u Nt exp [N t (4 u 4 — 4 u 6 + 0(u 8 ))] , 

with the fundamental coefficient of the character expan¬ 
sion uU3) = h{p)/h{p) = f (1 - £ + 0(/3 4 )). 

The contributions of the fermion determinant to the 
effective action S 2 are determined by hopping expanding 
the fermion determinant. Here we again follow the ex¬ 
ample of [22i and split the determinant in temporal and 
spatial hops, T = T + + T~ and S = S + + S~ in forward 
and backward direction, 

det [D\ = det[l -T+-T~ -S + - S~] 

= det[l - T] det[l - (1 - T)~\S+ + 5")] 

= det[D stat ] det[D kin ] . (7) 

The part containing only temporal hops, called the static 
part of the determinant, is easy to evaluate since we do 
not have to integrate over spatial links. By evaluating 
the spin and space determinant we get 

det[D stat ] = ]^[(1 + hLg + /i 2 ) 2 ( 1 + hLg + h 2 ) 2 , (8) 

X 

with 


with the fermion determinant det Df and the Wilson 
gauge action 

5 « = ^rE( trC/ P + trC/ p)- ( 2 ) 

c p 

The effective action is then defined by integration over 
the spatial links 

Z = j [dU 0 \ exp[—S e ff] , 

, N f 

exp[—5 e ff] = / [dUi\ exp[—^g] ]^[ det D f . (3) 

J /=! 


h = (2 n.e all ) Nt = and h = h{-fx) . (9) 

In the strong-coupling limit fj = 0 the connection be¬ 
tween the hopping parameter k and the constituent quark 
mass is then given by am q = — ln(2«) [30 . To evaluate 
the kinetic quark determinant we further split the deter¬ 
minant in into parts describing quarks moving in posi¬ 
tive Pi and negative M t spatial directions i. with P = 
E* Pi = (1 - T)- 1 ^ and M = E, Mi = (1 - T)~ 1 S~, 

det[D kin ] = det[l - (1 - T)~ 1 (S+ - S~)] 

= det[l - P-M] 

= exp[Tr log(l — P — M )] . (10) 


We can seperate the effective action 

S eS (W) = S l (W) + S 2 (W) , (4) 

where S\ contains all contributions from the Yang-Mills 
action potentially modified by contributions from non¬ 
winding fermion loops, while S 2 contains the winding 
fermion loops with gauge corrections. Both S\ and S 2 
now only depend on temporal Wilson lines or their traces, 
the Polyakov loops 


N t -1 

Lg = tr Wg = tr U 0 (x,t) . 
t =0 


Since the trace in |To| ) is also a trace in coordinate space, 
only closed loops contribute for which we need equal 
numbers of P’s and M’s. Up to order 0{k a ) this yields 
the following terms 


det [/-An]] — 


exp 


-Tr PM - Tr PPMM - -Tr PM PM 

2 


LV) 


= 1 - Tr PM - Tr PPMM - -Tr PM PM 

2 

+ ^(Tr PM) 2 + 0{k 6 ). 


(5) 


( 11 ) 
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When we evaluate the terms in we recover all the 
terms that are there in the case of SU( 3). However there 
are two additional terms that appear for SU( 2) in the 
expressions 

J dU y^tr PjMjPjMj , and 
\ J du Y, tr P Z Mi tr PiMi , (12) 

because of an additional non-vanishing group integral 

J dU UijUkl = tik^jl • (1^) 

We identify those terms with possible diquark contribu¬ 
tions in the two-color theory. With (13) the contributions 
from the diquark terms become: 

[ dUY^rPMAM, = 


16k ^ (](:‘i B xy ] , 




- / dU PjMjti PjMj = 


(14) 


32k 4 ^2 det i B x, y \ det [B x+i ^ y+i ] 


x,y 


where B XtV are parts of the static quark propagator, 


Atat — "4 + 7o B 


B X ,y ~ B^y + B x y 


(15) 


B+ = 

X ’V 2 1 + hW xv 


+ 2 Z 


1 W(t x , ty) 

1 + hW 


\0(t-y t x ) h6(t x ty)\ d xy ) 


1 hWi 

x ’ y ~ 2 1 + hWt 


1 = t x -t„ W^(t X , ty) 


+ -z** - v 


l + /iWt 


[6(t X - ty) - hO(ty - t X )\ 5 X 


and are defined in the same way as in [22] with z = 2 k e aM 
and z = 2Ke _a/i . Eqs. (141 then lead to an overall di¬ 


quark contribution of the form 

^diquark = 16 k ^ ) (let \B x _y. det \B x -\-i.y—/] (16) 


x,y 


= k 4 1V 2 J2 h 4 det 

x,i 

= « 4 a 2 £ 


1 


(1 + hWg)( 1 + hWs+i) 
h 4 


(1 + hLg + h 2 ){ 1 + hLg +i + h 2 ) 


strong-coupling description is still valid even at /? = 2.5. 
This is due to the fact that the effective coupling for the 
gauge sector A is negligible in the temperature regime we 
are interested in (T < 10 MeV). For the temperature to 
be this low we need a large temporal extend of the lat¬ 
tice, this leads to a negligible effective gauge coupling i.e. 
A(/3 = 2.5, N t = 200) ~ 1 • 10~ 15 . We therefore end up 
with a completely fermionic partition function. The only 
leftovers from the Yang-Mills part of the original theory 
come from gauge corrections to the effective fermion cou¬ 
plings, which are the same as in SU( 3) and amount to 
replacing h in (J9| by 


h = exp 


Nt ( afi T In 2k + 6k 2 


,n t 


1 — u 


■ (17) 


Now we can piece everything together to get the effective 
action for the cold and dense regime and one flavour: 


- S eS = ^2 l°s(l + hL i 
-2/^tr hW * 


-tr 


+ K 2 ) 2 

hWt 


x-\-i 


1 + hW$ 1 + hW^+i 


+ 2 


. k 4 1V 2 

N 2 

k 4 N 2 


■£* 


hW* 


rtr ■ 


hW; 


x-\-i 


N 2 


£t r 


(1 + hW s ) 2 (1 + hW s+i ) 2 
hW: 


tr -—tr J 




(1 + hW s ) 2 1 + hWg-i 1 + hW s ^ 


,K 4 iV 2 


hW* 


hW S -i hW 7 

rtr-——— tr 


■ 2 " ' r V tr 

N 2 (1 + hW s ) 2 “ 1 + hW S -i ” 1 + hW s+J 

x i^i3 




hW * 


hWg+i hWx-\-j 

rtr-——-tr - 


k 4 N 2 

N 2 “ (1 + hW s ) 2 “ 1 + hW s+l “ 1 + hW s+j 
h 4 

^ 5 (l + hL s + h 2 )(l + hL s+i + h 2 ) ’ (18) 


where the second fermion coupling h 2 is defined as 

,.N t 


h 2 = 


Nr. 


1 


1 — u 


(19) 


Note that the combined strong-coupling hopping expan¬ 
sion is defined as an expansion in K m u n , and given here 
to the order m + n = 4. Unlike the one-point fermion 
coupling h in (171, a partial resummation of h 2 is not 
possible at this order [22] • We furthermore dropped all 
terms proportional to h = h(~n) in Eq. (18) as well as 
terms that are subleading in N t . This implies that the 
form of the effective action in (18) is valid only for suffi¬ 
ciently large p. and N t . 


1. Strong Coupling Limit and Gauge Corrections 

When leaving the strong-coupling limit fj = 0 one 
might think that we should include Polyakov-loop inter¬ 
actions of the form @. However, it turns out that the 


2. Effective Action for Nf Flavors 

In the theory with Nf quark flavors we have to intro¬ 
duce fermion determiants for each of them in Eq. (|3|. 
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When they are degenerate this simply amounts to tak¬ 
ing the power Nf of the fermion determinant. As in the 
one-flavor case we split the fermion determinant into a 
static and a kinetic part and do the hopping expansion 
analogous to that in Eqs. ([7| and 


det [D] Nf =(det[Atat]det[Z? ki n]) iV/ (20) 

= det[D sta t] ^ exp \-N f tr PM - N f tr PPMM 

N f K i 

- tr PM PM + 0(k 6 ) 

= detfDstat]^ [l - N f tr PM - N f tr PPMM 

N f N 2 f „ 

- tr PM PM + (tr PM) 2 + 0(k 6 ) 


This leads to the effective action for Nf degenerate fla¬ 
vors: 


- S e ff = N f^2 log(l + hL s + h 


2\2 


( 21 ) 


— 2 Nfh ,2 tr 


hW* 


tr 


hW. 


x-\-i 


+ 2N 


x,i 

> ^N 2 
f N? 

k a N 2 


1 + hW s 1 + hW> 


E tr 


x-\-i 

hW* hW i 

t tr 


x-\-i 


(1 + hW s y (1 H- hW$+i) 2 


K v-^ 

+ JV /wfE 


tr 


hWg 


rtr 


hW*-i 


N 2 ^ (1 + hW s ) 2 ur 1 + hWg-i tr 1 + hW 3 -i 

X,l,J J 

k 4 7V 2 v-^ 

2> 




+ 2N 


hW 7 


tr 


hWm-i 


-tr 


hW. 


x+j 


N 2 (1 + hW s ) 2 1 + hWg-i 1 + hW, 




x+j 


Ar k a N 2 hWx 

+N f^ E tr •- 


hWx+i hWx-f-j 


(1 + hWg) 2 1 + hWg+i 1 + hW ; 


x+1 


+ ( 2iV / iV /) K ^ E (i + /ji 2 + /j2)(i + + h 2 ) ' 


B. Symmetries and Spectrum of two-color QCD 

One of the unique features of two-color QCD is the fact 
that the fundamental representation of SU(2) is pseudo- 
real, i.e. for vanishing chemical potential, quarks and 
anti-quarks belong to equivalent representations. This 
results in an enlarged flavor symmetry. In the contin¬ 
uum theory with Nf degenerate quark flavors the sym¬ 
metry group is SU(2Nf) which contains the usual chi¬ 
ral SU(Nf)L x SU(Nf)n x U(1)b symmetry as a sub¬ 
group, see e.g. m- For Nf = 2 continuuum flavors 
the enlarged flavor symmetry is 5Z7(4) and the chiral 
symmetry breaking pattern is 517(4) -+ Sp( 2) which 
is locally isomorphic to 50(6) —b 50(5). Dynamical 
chiral symmetry breaking then leads to five Goldstone 
bosons in the chiral limit: the three pions and a scalar 
diquark/anti-diquark pair. For a small finite quark mass 
m q these diquarks are thus pseudo-Goldstone bosons de¬ 


generate with the pions. As such they also are the light¬ 
est baryonic excitations in two-color QCD with bosonic 
baryons. The enlarged flavor symmetry is also broken 
by a finite baryon-chemical potential ps = 2/t > 0, 
as 517(4) -b SU(2) l x 517(2) fl x U( 1) B . With both, 
m q > 0 and ps > 0, one is left with the usual isospin 
and baryon number symmetries, SU(2)y x U(1)b- The 
two-color analogue of the nuclear-matter transition is 
Bose-Einstein condensation of the scalar diquarks, where 
a non-vanishing diquark condensate ( qq ) dynamically 
breaks the baryon number U(1)b- Because this tran¬ 
sition is of second order and hence continuous, for T = 0 
it must occur at hb = n^d, he. /t c = ntd /2 = rn^/2. 
In particular, there is no shift in /t c by binding energy 
as in the liquid-gas transition of nuclear matter in QCD 
(where 3p c = ms — e = 923 MeV for a nucleon mass 
ms = 939 MeV and a binding energy per nucleon of 
e = 16 MeV). 

From these general symmetry considerations, one 
would naively predict for Nf = 1 a chiral symme¬ 
try breaking pattern 517(2) —b 17(1 )b with a scalar 
diquark/anti-diquark pair as the two pseudo-Goldstone 
bosons. However, it turns out that one cannot construct 
a totally antisymmetric color-singlet pair of quarks in the 
scalar channel with only one flavor. There is therefore no 
diquark condesate for Nf = 1 either and baryon number 
remains unbroken, at least in the confined phase. 

Since the present approach is valid only for very heavy 
quarks, however, there is no approximate chiral symme¬ 
try in the first place F] Nevertheless, for Nf = 2 there is 
a scalar diquark which still remains exactly degenerate 
with the pion in the vacuum. It is therefore straightfor¬ 
ward to calculate its mass in the combined strong cou¬ 
pling and hopping expansion. The result is then the same 
as for the pion, 

amd = —21n(2tc) — 6k 2 — 24k 2 — 1 -b 6k 4 + • • • , (22) 

1 — u 

with corrections of the order n + m = 6. For com¬ 
parison, we can estimate the quark mass by assuming 
m q = /i — T In h as at leading order, but now with the 
strong-coupling and hopping expansion corrections in¬ 
cluded in h which are formally the same as in 517(3) and 
available to the order m + n = 7 [I5j. This corresponds 
to extracting the quark mass from the static propagator 
D~ t ^ t at this order and yields, 

am„ = — ln(2«) — 6k 2 —-- b 48k 4 u(1 — it) 

1 — it 

— 24 k 2 u (it 4 + k 2 u 2 — k 4 ) + ■ • • . (23) 

It amounts to a small binding energy for the scalar di¬ 
quark beyond leading order, in lattice units, 

2am q — amd = 6 k 2 y~— —6k 4 (1 — 8u(1 — it))+ • • • . (24) 


1 At least it means that we don’t have to worry about the poor 
chiral properties of Wilson fermions for which everything except 
the SU(2)y x U(1)b is broken explicitly by the Wilson term. 
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We reiterate that this scalar diquark does not exist in the 
one-flavor case. 


C. Leading-Order Mean-Field Density 


We can get some insight into the model by looking at 
the static part of the fermion determinant in the strong 
coupling limit A = /3 = 0. Then, the partition function 
factorizes, 


Z(t3 = 0) = (^J dW [1 + hLz + h 2 ] 2N f (25) 


x[l + hL s +h 2 ] 2N f) 


N“ 


For T —> 0 at a finite chemical Potential p, we have h —> 0. 
Within a mean-field description of the Polyakov loop this 
partition function then simplifies to 

Z = [1 + 2hL + h 2 ] 2NfN * , (26) 


where a factor of two was inserted because the mean- 
field Polyakov loop L here is normalized to assume values 
within [—1,1], The quark number density then follows as 


T d 

Vdjt 1 2 


a 3 n = 4V/- 


1 + Le t 


1 + 2 Le t + i 


-M) 


(27) 


The same expression can be derived from a Polyakov- 
Quark-Meson-Diquark model 1291 for very heavy quarks. 


For L = 1 it simply reduces to the zero-momentum oc¬ 
cupation number of the Fermi-Dirac distribution. In the 
deconfined phase it thus describes a free gas of heavy 
quarks. We will see below that for small but finite 
Polyakov-loop expectation values L > 0, and small T, 
Eq. (271 due to imperfect confinement behaves as a quark 
gas, suppressed by L 1 up to some critical chemical po¬ 
tential n c (T) where it starts to reflect the behavior of a 
diquark gas with m<i = 2 m q at this leading order. Explic¬ 
itly, consider /i < m q and hence x = exp {{m q — u)/T} > 
1. We can then distinguish two regimes: 



FIG. 1. Comparison of quark densities evaluated at the orders 
k 2 and fid, both for k = 0.00802, /3 = 2.5, fx = uid /2 and N t 
ranging between 100 and 2400, including gauge corrections. 


seen from Equation (18) the relevant expansion parame¬ 
ter for the effective theory is n 2 N t /N c . Because we are 
interested in very low temperatures T = l/aN t , we need 
lattices of large extend N t in the temporal direction, es¬ 
pecially as we go to smaller lattice spacings a. Therefore 
our hopping parameter n needs to be sufficiently small. 
This implies that one can reach smaller quark masses at 
higher temperatures and vice versa. To check the conver¬ 
gence of the hopping expansion we compare expectation 
values for the densities a 3 n from simulations including 
corrections up to 0(n 2 ) and 0(k 4 ) for different values 




of the expansion parameter Figure [1 shows a plot 

of the two densities for /i = Wd/2 at /3 = 2.5, including 
gauge corrections. The two agree reasonably well for all 
values of n 2 N t /2 up to slightly above 0.02. 


Scale Setting and Units 


Assuming that quarks as heavy as those used here only 
have a negligible influence on the running of the coupling, 
we use the non-perturbative /3-function from |32j 


1 

ay/a 



(29) 


a 3 n 

4 Nf 


e (nB-m d )/T ) Xjc<1, 

L e (n-m q )/T 5 Xar» 1 . 


(28) 


with the parameters 

d= 1.98(1), 6 = 0.305(6). 


(30) 


The transition occurs at around fi c ~ m q + T In L . 

III. RESULTS 

In order to assess the range of validity of the effective 
theory we need to test the convergence properties of the 
hopping series for various parameter choices. As can be 


In lack of phenomenological input we somewhat arbitrar¬ 
ily use the typical yfa = 440 MeV also for the string 
tension of two-color QCD. With Eqs. (29) and (22) we 
are then able to assign physical scales to our systems. In 
all our simulations the diquark mass from Eq. (22) is ad¬ 
justed in this way to vrid = 20 GeV, and the temperatures 
range between T = 3.454 MeV and 9 MeV. On our finest 
lattice with /3 = 2.5, corresponding to a = 0.0810 fm, this 
amounts to k = 0.00802123 and N t values between 269 
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9.9 9.92 9.94 9.96 9.98 10 10.02 


l_i in GeV 

FIG. 2. Quark density a 3 n in lattice units and Polyakov-loop 
expecation value (|L|), both over /x, with simulation parame¬ 
ters f3 = 2.5, k = 0.00802, leading to m,d = 20 GeV, N s = 16, 
Nt = 484, corresponding to T = 5 MeV, and Nf = 1. 



FIG. 3. The pseudo-critical line for the deconfinement 
transition from simulations with the same parameters as 
in Fig. [2] except for different Nt, corresponding to T = 
9, 8,... 3.454 MeV, and a linear extrapolation to T = 0. 


and 700 such that n 2 N t /2 lies within 0.0087 and 0.0225, 
and hence within the range of validity of the hopping 
expansion, cf. Fig. [l] The parameters for the coarser 
lattices with f3 values down to 2.4 lead to even smaller 
values for the expansion parameter. 


A. Results for Nf = 1 

In this subsection we first present our numerical results 
for the Polyakov loop and the quark density in the effec¬ 
tive theory at low temperatures with one quark flavor, 
N f = !• _ 

As an indication of deconfinement at high density and 
low temperature we plot the the Polyakov-loop expec¬ 
tation value in Fig. [2j where (|L|) stands for the usual 
expectation value of the modulus of the volume averaged 

L = — L$ . (31) 

X 

Because of the presence of dynamical quarks, even for 
vanishing net-baryon density, the Polyakov loop will have 
a small but nonzero expectation value (|L|) > 0. For the 
parameters of Fig. |2] its /x = 0 value is (|L|) = 0.012. For 
aspect ratios as the one considered here, with N t /N s « 30 
in this example, this value is determined by the finite 
spatial volume. It is therefore basically temperature in¬ 
dependent. It furthermore also remains constant in /i 
until just below the onset of the density near m ( //2 be¬ 
cause the temperature of T = 5 MeV here is so low that 
no baryonic degrees of freedom are being excited as long 
as the baryon chemical potential 2/x stays well below the 
gap in the baryon spectrum. From Fig. [2] it appears how¬ 
ever that (\L\) starts to rise from its /x = 0 value before 
the onset of the density so that we can not distinguish 
baryon density from quark density here. As we will dis¬ 
cuss in the next subsection, however, there is a very small 


diquark contribution in the density above /x c ss 9.96 GeV 
at this temperature which is not visible on the linear 
scale of Fig. [5] but where the Polyakov loop still has its 
zero-density expectation value. We will provide some ev¬ 
idence that this diquark-density onset might stay below 
the deconfinement transition when we extrapolate both 
to T = 0 for N f =2 with scalar diquark below. 

At larger values of the chemical potential the quark 
density saturates at a 3 n = 2 N c Nf, the maximum num¬ 
ber of quarks per site due to the Pauli principle, as in the 
effective theory for heavy quarks in SU( 3) [ST]. This be¬ 
havior which is a lattice artifact has previously also been 
observed in finite density simulations of two-color |26] 
and G 2 -QCD [SB:- This saturation leads to an effective 
quenching of the quarks and hence the Polyakov loop 
decreases again as it is approached. We also analyzed 
the Polyakov-loop susceptibility and found no increase of 
its rather broad maximum with the lattice volume hence 
indicating a smooth cross-over behavior rather than a de- 
confinement phase transition in the infinite volume limit. 
The pseudo-critical chemical potential /x pc from the in¬ 
flection point of (|L|) along the /x axis is shown in Fig. [3] 
It coincides with the point where (|L|) reaches half its 
maximum value. By determining this /Lt pc for different 
temperatures we obtain a pseudo-critical line for the de- 
confiment transition at low temperature which can be ex¬ 
trapolated to T = 0 as shown in Fig. [3] We can see that 
the pseudo-critical line terminates at /x pc slightly above 
10 GeV which corresponds to half the scalar-diquark 
mass from Eq. (22) with the parameters used here. 

We conclude this subsection by discussing in some 
more detail the quark-number density n defined by 


n =vb' nZ - (32) 

As mentioned above, Fig. [2] shows the strong increase in 
the density in lattice units a 3 n at a value of the quark 
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FIG. 4. Comparison of the quark densi ty a t order k 4 from 
Fig .[ijto the leading-order formula in Eq. ( |27| with L replaced 
by the corresponding data for (\L\)/2 from Fig. [2] T = 5 MeV. 


FIG. 5. Logarithmic plot of the density in lattice units a 3 n 
and the Polyakov loop (|L|) as a function of p in physical 
units at T = 5 MeV with the same parameters as in Fig. [2] 


chemical potential around /z = m d / 2, from where on it 
rapidly grows to its saturation value with each lattice site 
fully occupied by 2N c Nf = 4 quarks. This transition is 
described qualitatively well by the leading-order mean- 
field formula in Eq. (271. This can be seen in Fig. [4]where 
we compare the data for the quark density at order n 4 
from Fig. [2]to the leading order form in Eq. ( p7| with the 
corresponding am q = — ln(2«) for the strong-coupling 
limit, and with L replaced by (|L|)/2, i.e., using the /z- 
dependent data for the Polyakov-loop expectation value 
of Fig. [2] in the mean-field approximation. 

To resolve the differences we need to have a closer look 
at the behavior of the chemical-potential and tempera¬ 
ture dependence of the quark density, especially in the 
region where L x ~ 1 with x = exp{ ( m q — /z) /T } as de¬ 
fined in Sec. |II C| above. Fig. [5] shows a logarithmic plot 
of the density and the Polyakov loop of Fig. [2j 

We observe two different regimes of exponential in¬ 
crease before the density approaches its saturation value. 
They are separated by a kink in the logarithmic plot, 
here at /z c « 9.96 GeV, where the Polyakov-loop still is 
constant at its /z = 0 expectation value, (|L|)/2 ~ 0.006. 


The two regimes correspond to the two limits in (28) 


of the leading-order mean-held density, Eq. (27). This is 
demonstrated in Fig. [6] where we compare the density to 
two corresponding fits: 

When we fit the data in the region of the second ex¬ 
ponential increase, for fi values between 9.96 GeV and 
9.99 GeV, to 


a 3 n = 4Nf exp{(2/z — mat)/T} 


(33) 


with a single parameter mat we obtain for Nf = 1 , 

mat = 20.0045(5) GeV. (34) 

For comparison, with the same lattice parameters the 


quark mass from Eq. (|23j) becomes m q = 10.0024 GeV 


at the leading order n + m = 0, niq ’ = 10.0014 GeV at 

(7) 

the order n+m = 4, and m q = 10.0013 GeV at the order 
n + m = 7. Therefore, the fit parameter mat is consistent 
with 2mj°' but slightly larger than 2 m q = 20.0028(2) at 
the same order n + m = 4 (with an error of the size 
of the higher-order corrections up to n + m = 7 as given 
explicitly in Eq. (23)). It is larger than the corresponding 


scalar diquark mass, m d — 20 GeV from Eq. (22), which 


might simply reflect the fact this scalar diquark does not 
exist for Nf = 1. A daring interpretation would be that 
we see a heavier diquark mass here, such as that of an 
axial-vector diquark which can exist also for Nf = 1 . 

To test this, we have done the same analysis with the 
same parameters also for two flavors (see the next sub¬ 
section). The same fit to the form in (33) then yields for 

Nf = 2, 


mat = 19.9986(10) GeV, 


(35) 


which is now indeed very close to the scalar diquark mass 
md and significantly smaller than 2 m q . In order to quan¬ 
titatively describe this regime of exponential increase we 
therefore have to replace 2 m q by in the leading-order 
mean-held formula for the density in Eq. (27). This in¬ 


dicates that matter on this side of the kink consists of 
diquark excitations. 

The first exponential increase, for the /z values below 
9.96 GeV, is described by 


a 3 n = ANf L exp{(/u — m q )/T} . 


(36) 


In this case we use m q = 10.0014 GeV for the quark 
mass at this order and fit the data via L as the free 
parameter. For the one-flavor data of Fig. [6] this leads 
to L of the order of 10 —4 , however, with a very large 
uncertainty. It determines the precise value of the onset 
He of the diquark density, by the intersection point of the 
two different exponential fits (33) and ((361), as 


L = exp{(/z c -m d + m q )/T} . 


(37) 




















^ in GeV 

FIG. 6. Logarithmic plot of the density in lattice units a 3 n 
at T = 5 MeV compared to one-parameter fits using L for 
fj, < 9.96 GeV and rrid for p > 9.96 GeV, see text. 



FIG. 7. Comparison of the measured density with the leading- 
order mean-held density with L = (Lg)/2 at mj = 20 GeV, 
j3 = 2.5, k = 0.00802, and different temperatures. 


With nid = mat from (34) for the Nf = 1 data in Fig. [6j 
for example, this leads to values between L = 8 • 10 -5 for 
H c = 9.956 GeV and L = 1.2-10" 4 for /. t c = 9.958 GeV. In 
any case, it is much smaller than the zero-density value of 
(\L\)/2 « 0.006. Instead we observe that it is more con¬ 
sistent with the expectation value of the local Polyakov 
loop Lg. Its expectation value is extracted from the per- 
site probability distribution P(Lg) which we otain by 
histograming the local Polyakov-loop variable Lg as in 
[32] . At T = 5 MeV, with /3 = 2.5 and k = 0.00802, we 
obtain for this observable a zero-density value of about 
(Lg) ~ 10~ 4 instead of (|L|) ss 0.012 for the modulus 
of the volume-averaged Polyakov loop. This suggests 
that one should use the local Polyakov-loop expectation 
value (Lg) in mean-field approximations as Eq ( |27[ ) . Like 
(|L|) it is independent of the chemical potential below 
the deconfinement crossover at n pc . And as soon as the 
Polyakov-loop starts to rise from its constant zero-density 
expectation value the two agree well within the errors. It 
is only the residual small value at vanishing net-baryon 
density due to imperfect confinement in a finite volume 
in which the two differ simply because it takes longer for 
the modulus to vanish than the local Polyakov loop in 
the infinite-volume limit for /r < fi pc . This difference is 
only relevant at densities in lattice units below 10 -4 and 
hence not visible on the linear scale of Fig. [4] above. 

In fact, using the local Polyakov-loop expectation value 
(Lg) for various temperatures in the leading-order mean- 
field formula for the quark density, Eq. (l27|), describes 


the data especially also in the low-density region around 
the diquark-density onset at /j, c very well as can be seen 
in Fig. [?] These are not fits. We simply use Eq. (27) 


with L = (Lg) /2 here to describe the quark density over 
the whole range of temperatures we have investigated. It 
describes the imperfect statistical confinement of quarks 
for p, below /z c = — m q + Tln((Lg) /2), an ensemble 

of diquarks above /z c , and quark matter with lattice sat¬ 
uration for /r larger than /i pc where 2 L = (Lg) = (|L|) as 



1.2 


0.2 


H in GeV 


FIG. 8. Quark density in lattice units and Polyakov-loop 
expectation value (|L|), both over fj,, with the same simulation 
parameters as in Fig. [2j but for Nf = 2 here. 


in Fig. [4j all at the same time. 


B. Results for Nf = 2 

In this subsection we discuss the results from the ef¬ 
fective theory with Nf = 2 degenerate quark flavors in 
somewhat more detail. In particular, we describe how 
well the onset of the diquark density at fL = fi c agrees 
with the scalar diquark mass md from Eq. (22) which we 
know that it exists for Nf = 2. 

As a first example in Fig. IS] we present the Nf = 2 
results for the quark density and the Polyakov loop on our 
finest lattice, with /3 = 2.5, k = 0.00802123, N s = 16 and 
N t = 484 corresponding to md = 20 GeV and T = 5 MeV 
at a lattice spacing of a = 0.081 fm as in the previous 
subsection for the one-flavor case. Both observables show 
the same qualitative behavior as for Nf = 1 before. 
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FIG. 9. Comparison of the Polyakov loop (|L|) for Nf = 1 
and 2 at T = 5 MeV with parameters as in Figs. [2] and [8] 



FIG. 11. Section of the Nf = 2 phase diagram with a linear 
extrapolation of the diquark-density onset to T = 0 at /3 = 2.5 
or a = 0.081 fm with N„ = 16 and N s = 48. 



FIG. 10. Difference between the normalized quark number 
densities for Nf = 2 and A r / = latT = 5 MeV. 


The most obvious differences between Nf = 1 and 2 
are the different saturation densities given by a 3 n sa t = 
2 N c Nf and the maximum value of the Polyakov loop 
(|I/|) which is a bit higher for Nf = 2. The direct com¬ 
parison of the Polyakov-loop expectation values in Fig. [9] 
shows that the deconfinement crossover tends to start at 
somewhat smaller values of fi for Nf = 2, and it gets 
quenched later, likewise. 

The difference between the normalized quark-number 
densities n/n sa t for Nf = 2 and Nf = 1 is shown in 
Fig. IT We can see a deviation around fi = 10 GeV. This 
is in line with the observation that the deconfinement 
transition happens earlier for Nf = 2 as well. 

As described for Nf = 1 above, we follow the same 
procedure with Nf = 2 for various temperatures from 
9 MeV down to 3.454 MeV. That is, for each tempera¬ 
ture we determine the intersection point of the two expo¬ 
nential regimes in the quark density (for Nf = 2 their fi- 
values are consistently about 1 — 2 MeV lower than those 
for Nf = 1). Since the mass-parameter in the second 


exponential for Nf = 2, see Eq. (35), agrees well with 
the scalar diquark mass, md = 20 GeV from Eq. (221, 
we now take the intersection of the lines in logarithmic 
plots analogous to Fig.[6]as the onset of baryonic diquark 
matter and extrapolate the corresponding onset chemi¬ 
cal potentials fi c to T = 0. The result for the N s = 16 
lattice is shown in Fig. 0 Using a linear extrapolation 
as in the figure, which is consistent with a temperature 
independent L ss {L)/2, the result for the T = 0 diquark 
onset on the N s = 16 lattice is fi c = 9.9996(22) GeV and 
hence includes mj /2 = 10 GeV within the error. Larger 
lattices lead to smaller values of L and hence a smaller 
slope In L in the linear extrapolation, but the extrapo¬ 
lated fi c remains the same. For comparison, the same 
analysis was also done on a N s = 48 lattice with the 
result that fi c = 9.9998(9) GeV as also shown in Fig. [6} 
In order to test the scaling of this onset we have per¬ 
formed the same analysis also for 7 different lattice cou¬ 
plings /3 between 2.4 and 2.5, corresponding to lattice 
spacings between a = 0.1124 fm and 0.0810 fm with n 
values adjusted so that nid from Eq. (22) remains fixed 
at 20 GeV as before. Again, for each /3 we extract the 
corresponding intersection points of the two exponential 
regimes in the quark density at the same 7 temperatures 
between 9 MeV and 3.454 MeV. The extrapolated Nf = 2 
results for the zero-temperature diquark-density onsets 


from these intersection points are collected in Fig. 12 


Within the errors, these extrapolated values for fi c ba¬ 
sically all agree. Assuming that fi c is thus independent 
of the lattice spacing in this parameter regime we simply 
use their average as our final overall estimate of 


fi c = 10.0001(3) (38) 

from the data in Fig. [12] as indicated by the horizontal 
line with the gray error band. This overall estimate thus 
confirms that fi c = md/2 with rather high precision. 

This agrees with the corresponding onset of isospin 
density at m ^ r /2 in the effective theory for heavy quarks 
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FIG. 12. Critical chemical potentials fi c for the diquark- 
density onset at T = 0 from linear extrapolations for different 
lattice spacings between a = 0.0810 fm and 0.1124 fm corre¬ 
sponding to lattice couplings between /? = 2.5 and 2.4. 


in QCD [ 22 ] j and it shows that there is no “Silver Blaze” 
problem [33] in the effective lattice theory for two-color 
QCD with heavy quarks either. Unlike the nuclear- 
matter transition at 3/t c = tob — e, with some evidence 
of a finite binding energy per nucleon e in the effective 
theory for heavy quarks in QCD [22] . there is certainly 
no such evidence of a shift of the onset /i c by a non-zero 
e here. This is consistent with the generally expected 
difference between a first-order liquid-gas transition in 
QCD and Bose-Einstein condensation of diquarks in a 
second-order transition in two-color QCD. At the same 
time, however, the diquark densities obtained here are 
far from reflecting any sign of Bose-Einstein condensa¬ 
tion. Quarks and cliquarks are way too heavy to inter¬ 
pret the latter as deeply bound dimers. With our pa¬ 
rameters from Eq. (241 the 20 GeV diquarks are only 
bound by about 2.8 MeV. If it wasn’t for confinement, 
the transition temperature of the diquark-condensation 
phase by pair breaking should roughly be of the same or¬ 
der. Probably only because of the statistical confinement 
of the quarks in the first place, at the available temper¬ 
atures above 4 MeV, all we can observe here therefore 
is an essentially free heavy-diquark gas behavior in the 
small window between fi = fJ. c {T) and the beginning de¬ 
confinement crossover followed by lattice saturation. Un¬ 
fortunately, the region where one might find a superfluid 
diquark-condensation phase thus is currently still beyond 
reach within the convergence region of the hopping se¬ 
ries. Nevertheless, we can attempt to give a very rough 
first estimate of a region where such a diquark superfluid 
might be found, if we were able to further reduce the 
temperature, as follows: 


Since the pseudo-critical chemical potentials for decon¬ 
finement at the available temperatures are also all below 
md/ 2, we compare their zero-temperature extrapolation 
to that of the diquark-density onset in Fig. [13] The dif¬ 
ference between the so extrapolated /r c = Wd/2 and n pc is 



9.985 9.99 9.995 10 

|j in GeV 

FIG. 13. Extrapolation of the deconfinement transition com¬ 
pared to the extrapolated diquark-density onset from Fig. m 
The region where one might hope to find a superfluid diquark- 
condensation phase is marked by the shaded red triangle. 
Half the diquark mass, m,j/2 = 10 GeV, and the quark mass 
m q = 10.0014 GeV are indicated by vertical lines. 


small but significant. As seen in the figure, the deconfine¬ 
ment crossover then hits the T = 0 axis of the phase di¬ 
agram just above m q = 10.0014 GeV. Therefore, a small 
window for a potential superfluid diquark-condensation 
phase at sufficiently low temperatures remains. The re¬ 
gion where this might occur is indicated by the shaded 
red triangle in Fig. [T3] This region starts at a chemi¬ 
cal potential slightly below /i = md /2 = 10 GeV, i.e. at 
the lower limit given by the extrapolation error of the 
diquark-density onset for the f3 = 2.5 data used here. We 
use this lower limit instead of to^/ 2 because there are also 
some truncation errors in our equations for the diquark 
mass m.d , Eq. (22), and the effective fermion couplings h 
and ft -2 in Eqs. (17) and (19). The deconfinement tran¬ 
sition temperature at this lower limit is then of roughly 
the same order as the diquark-binding energy and hence 
of the naive estimate of the transition temperature of a 
possible diquark-condensation phase. 


IV. SUMMARY AND OUTLOOK 


In this paper we have presented results for the bary- 
onic diquark-matter onset in two-color QCD with heavy 
quarks. The results were obtained by numerical lat¬ 
tice simulations of a three-dimensional effective theory 
of Polyakov-loop variables analogous to that previously 
used for QCD [ 22 ] • The effective theory is thereby de¬ 
rived from the full theory by a combined strong-coupling 
and hopping expansion, in this case up to the combined 
order in ( u n n m ) with n + m = 4. This effective theory 
was applied to the cold and dense regime of the phase di¬ 
agram of two-color QCD with both TV/ = 1 and TV/ = 2. 
We have mapped the pseudo-critical line of the decon¬ 
finement crossover at small temperatures as indicated 
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by the Polyakov-loop expectation value. With the scale 
set by the vacuum string tension as in QCD, we could 
reach temperatures ranging between 3.5 MeV and 9 MeV. 
At these low temperatures the Polyakov-loop expecta¬ 
tion value (L) stays constant at its small but finite zero 
net-baryon density value for chemcial potentials up to 
the deconfinement crossover at around a pseudo-critical 
/r P c- Within this constant regime we could clearly iden¬ 
tify a kink in the logarithm of the quark-number density 
where its finite-temperature behavior changes between 
that of a heavy-quark gas with imperfect statistical con¬ 
finement in a finite volume and an essentially free gas 
of heavy diquarks with mass rrid in a small window of 
chemical potentials // between this diquark-density on¬ 
set at n c ~ m<j/2 + Tln(L) (here with (L) £ [—1,1]) and 
/i pc . For fi > n pc the density describes quark matter that 
approaches its saturation value as the lattice gets filled 
with the maximum number of quarks per site allowed by 
the Pauli principle, and as expected from previous finite- 
density studies [2U1231HB] • 

This is all described quite well already by the leading- 
order analytic formula for the density together with a 
mean-field description for the Polyakov loop provided we 
make the following two adjustments: 

The constant Polyakov-loop expectation value (L) 
needed to describe the heavy-quark gas with imperfect 
statistical confinement for /i < fi c , below the diquark 
onset, is only consistent with the very small but finite 
expectation value of the local Polyakov-loop variable Lg 
as obtained from its unquenched probablility distribution 
which is slightly distorted by the large but finite quark 
mass mn q . This local Polyakov-loop expectation value is 
about two orders of magnitude smaller in this constant 
regime than the expectation value of the modulus of its 
volume average (\L\) usually used as the order parameter. 

In the heavy-diquark gas regime we have to replace 
2 m q by as computed analytically in the combined 
strong-coupling and hopping expansion. Especially for 
Nf = 2 the relevant mass rrid here is that of a scalar 
diquark degenerate with the pion. With our parameters 
we only have 2 m q — m,d = 2.8 MeV as the binding energy 
of this scalar diquark but that makes a significant differ¬ 
ence in the exponential increase of the diquark density 
for /j, c < /j, < n pc . 

As soon as the Polyakov-loop starts to rise from its 
zero-density value, {Lg) and (|L|) approach one another 
very rapidly. For // > \i pc we have {Lg) = (|L|) and their 
common //-dependent value in the leading-order density 
describes the approach towards lattice saturation quali¬ 
tatively quite well. 

The leading-order formula for the density can also be 
calculated analytically without the mean-field prescrip¬ 
tion for the Polyakov loop, of course. Then, however, it 
has only a single exponential increase and there is no di¬ 
quark onset just as there is no difference between nid and 
2 m q . Furthermore, one can also calculate the Polyakov- 
loop variable L in the mean-field formula from the exact 
leading order density for comparison. Its maximum at 


H = m q = — ln(2/c) for Nf = 1 then only results to be 
L max = -y/17/16 — 1 ~ 0.031, for example, which is more 
than an order of magnitude smaller than that of the mea¬ 
sured Polyakov-loop expectation value. 

On the other hand, with the adjustments mentioned 
above the leading-order mean-field formula is flexible 
enough to describe the G(« 4 )-data. The agreement be¬ 
tween the measured quark-number density and the mean- 
field treatment of the Polyakov loop in the leading-order 
form might in fact be somewhat surprising, because the 
measured Polyakov-loop distribution is not at all sharply 
peaked around the mean-field value. It is a rather 
broad distribution suggesting non-negligible fluctuations 
around the mean value. A similar effect for an effective 
Polyakov-loop model with dynamical fermions in SU{3) 
was found in |34j . 

The relatively detailed understanding of the quark 
number density allows us to interpret the most inter¬ 
esting baryonic regime fi c < /i < fi pc as an essentially 
free heavy-diquark gas with some confidence. There is 
no direct evidence of those baryonic diquarks forming a 
Bose-Einstein condensate (BEC), however. The temper¬ 
atures from 3.5 MeV upwards are still too high for our 
20 GeV diquarks bound by only about 2.8 MeV. Our rel¬ 
atively stable T —> 0 extrapolations of diquark-density 
onset and the deconfinement crossover nevertheless al¬ 
lowed to identify a small region between the two where 
one might hope to find signs of a BEC in the future. 

The fact that the T = 0 extrapolation of the diquark- 
density onset n c is scale independent and from our over¬ 
all average, covering various different lattice couplings 
between j3 = 2.4 and 2.5, agrees with mj/2 within a rel¬ 
ative accuracy of 2 • 10~ 5 at least is a strong indication 
that the zero-temperature diquark matter forming above 
fi c = irid /2 in two-color QCD does not involve binding 
energy which would shift the onset to values smaller than 
the lightest baryon mass as in QCD. In the future one 
might think about including a diquark source term and 
analyzing also a possible diquark condensate. In order to 
assess how much we can stretch the effective theory to¬ 
wards lighter masses and/or lower temperatures, which 
both increases the expansion parameter of the effective 
theory, it might be instructive to compare our results 
to full two-color QCD simulations with heavy quarks at 
finite //. Analogous work for G 2 -QCD with fermionic 
baryons and evidence of a first-order liquid-gas transi¬ 
tion of G 2 -nuclear matter [35] is currently in progress. 
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